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Abstract 

Atom interferometers using Bose-Einstein condensate that is confined in a waveguide and manip- 
ulated by optical pulses have been limited by their short coherence times. We present a theoretical 
model that offers a physically simple explanation for the loss of contrast and propose the method 
for increasing the fringe contrast by recombining the atoms at a different time. A simple, quantita- 
tively accurate, analytical expression for the optimized recombination time is presented and used 
to place limits on the physical parameters for which the contrast may be recovered. 
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I. INTRODUCTION 



The realization of a practical chip-based atom interferometer, using a Bose-Einstein con- 
densate, (BEC) would revolutionize internal navigation systems, precision measurements, 
and, perhaps, quantum information technology. The production and manipulation of a con- 
densate using an atom chip is appealing for many reasons. The trapping potentials, for the 
same electrical current, have a much higher frequency then their larger scale counterparts 
[1]. As a result, these devices require less power. An atomic cloud confined in a trap with 
a large frequency can be cooled to ultra-cold temperatures more rapidly [2], making the 
realization of a quasi-continuous source of condensed atoms more feasible. Chips also enable 
integration of several relatively simple components into devices with complex functionality 
[3, 4]. Finally, the entire device, including the vacuum chamber, can be miniaturized [5], 
making chip-based devices the most likely path to practical applications. 

The first experimentally realized BEC-based interferometers used a double- well potential 
to manipulate the atoms [6, 7]. In this type of interferometer, atoms are initially prepared 
in a single potential well that confines the atoms in all three directions. Then the well is 
continuously transformed into a symmetric double well, which must be done adiabatically 
slow to avoid collective excitations. After the sphtting, the wells are physically separated 
and, due to differences in the local environment, a quantum phase shift may develop between 
the atoms in each well. Nonlinearities, caused by atom-atom interactions, usually cause 
problems at the recombination stage in these types of interferometers. As a result, the 
atoms are recombined by suddenly switching off the trap allowing the two condensates to 
ballistically expand, overlap, and interfere. The phase shift can be measured by fitting the 
atomic density to a cos function. 

An alternative method for realizing a BEC-based atom interferometer uses a standing 
light wave created by two counter-propagating laser pulses to manipulate the condensate 
[8-10]. The BEC is loaded into a waveguide providing confinement along two directions 
but not in the third one. The standing hght wave splits the BEC into two harmonics 
moving in opposite directions. A second laser pulse reflects the harmonics by reversing their 
momentum. Finally, a third pulse recombines them, thus completing the interferometer 
sequence. Since the atoms are recombined in the guide, the phase shift between the two 
arms of the interferometer can be measured by counting the atoms in the zero momentum 
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state at the end of the interferometer cycle. 

In the first experimental realization of this type of interferometer [8] , the coherence time 
was about 10 ms. The loss of coherence was theoretically explained as being due to dis- 
tortion of the phase across each atomic cloud [11], which is caused by both the atom-atom 
interactions and the residual external parabolic potential in the axial direction. 

A similar interferometer with a coherence time of 44 ms has been realized in Ref. [9]. 
The coherence time of this interferometer was increased in two ways. First, the confining 
waveguide was created using macroscopic conductors, instead of microscopic conductors 
in a chip-based device. This larger scale device had looser confinement thus reducing the 
effect of the atom-atom interactions. Additionally, the waveguide was farther from the 
conductors minimizing small perturbations due to defects in the conductors or instabilities 
in the current. Second, the interferometer sequence used two reflections, such that each half 
of the condensate went through the same path twice in opposite directions. This method 
greatly reduces the phase distortion caused by the residual axial potential. 

A third interferometer was recently realized in [12]. This chip-based device has a coher- 
ence time of about 15 ms and uses a sightly different splitting technique. One of the lasers 
is frequency-shifted with respect to the other laser, resulting in a traveling wave optical po- 
tential. Using a 7r/2 pulse, the condensate is split into two harmonics one of which remains 
at rest and the other propagates with the momentum 2hki. A tt pulse acts as a mirror, and 
finally a second 7r/2 recombines the condensate. As a result of this different splitting mecha- 
nism, the splitting and recombination does not take place in the same physical location and 
the device is an atom Mach-Zehnder interferometer. By varying the radial confinement of 
the waveguide, the researchers were able to change the effective strength of the atom-atom 
interactions. They used this to test the theoretical model of Ref. [11] and demonstrated its 
validity. 

In this paper, we discuss possible ways to increase the coherence time of an atom inter- 
ferometer using Bragg diffraction to manipulate the atoms. Simple analytical expressions 
for the dynamics of the interference sequence are derived in the framework of the hydrody- 
namic approximation. Their validity is confirmed by their comparison with direct numerical 
solution of the Gross-Pitaevskii equation. Wc demonstrate that the coherence time can be 
increased by changing the recombination time and present simple analytic expressions for 
the optimized recombination time and the contrast. 
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In the rest of the section, we present a simple physical explanation for the loss of contrast 
in a BEC Michelson interferometer and discuss possible ways of restoring it. 

The interferometric cycle of duration T starts by illuminating the motionless BEC cloud 
with a splitting pulse from a pair of counter- propagating laser beams. This pulse acts like 
a diffraction grating splitting the cloud into two harmonics and ip^. The atoms diffracted 
into the +1 order absorb a photon from a laser beam with the momentum hki and re-emit 
it into the beam with the momentum —hki acquiring the net momentum 2hki. The cloud 
■0+ starts moving with the velocity Vq — 2hki/M, where ki is the wavenumber of the laser 
beams and M is the atomic mass. Similarly, the cloud ip- starts moving with the velocity 
—Vq- The two harmonics are allowed to propagate for the time T/2 and are illuminated by 
a reflection optical pulse. The atoms in the harmonics ijj^ change their velocity by —2vq 
and those in the harmonics -0- by 2vo. The harmonics propagate back for time T/2 and are 
subject to the action of the recombination optical pulse. After the recombination, the atoms 
in general populate all three harmonics ■0o and ip±. The degree of population depends on 
the relative phase between the harmonics ip± acquired during the interferometric cycle and 
can be used to deduce this phase. In particular, the wave function of the zero-momentum 
harmonics i/jq after the recombination is equal to 

iPo^^{ip+ + ip-), (1) 

where xjj± are the wave functions of the ±1 harmonics immediately before the recombination. 

Because of the nonlinearity and/or the external potential, the harmonics il)± do not travel 
with the velocities ±Vq during the cycle. First, the cloud "climbing up" the external potential 
slows down and the one moving "downhill" speeds up. Second, because of the nonlinearity, 
the speeds of the two clouds after their separation will be slightly larger than vq if the 
nonlinearity is repulsive and slightly less than vq if it is attractive. For definiteness, we 
shall discuss the influence of the repulsive nonlinearity assuming that the external potential 
is zero. An ideal operation of the interferometer in this case corresponds to all the atoms 
populating zero-momentum harmonics ipQ after the recombination, i.e., to Nq = Nfot- 

Because of the atom-atom interaction, the clouds iJj± exert a repulsive force on each other 
during the time they overlap. This force accelerates each cloud so that after the separation 
pulse the ip± harmonics propagate with velocities ±(i>o + Sv), where Sv > 0. The reflection 
pulses impart the momenta ^4:hki to the clouds transforming their roles: ip± — '0=f- After 
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the reflection, the ±1 harmonics propagate with the velocities ±(i>o — Sv). Harmonics' 
deceleration due to mutual repulsion during their overlap decreases the velocity of each 
harmonics by an additional 5v so immediately before the recombination the harmonics' 
velocities are ±(t>o — 26v). 

The nonzero value of Sv results in two consequences. First, since the clouds' speeds 
after the reflection pulse are smaller than before the pulse, the ±1 harmonics at the nom- 
inal recombination time still do not overlap each other completely. This effect is typically 
not very signiflcant. Much more important is the fact that the returning harmonics have 
momenta that are not equal to ±2hki and therefore can not be compensated by the recom- 
bination pulse. As a result, the wave function of the zero-momentum harmonics ipQ after 
the recombination can be written as 



where Ak — 2M5v/h and n{x) are the density proflles of the harmonics; their possible 
incomplete overlap has been neglected. The population of the zero-momentum harmonics 
is obtained by the spatial integration of {ipol"^- For AkR <C 1, where R is the characteristic 
size of the clouds, all the atoms after the recombination are indeed in the zeroth harmonics, 
i.e. A^o — ^tot- In the opposite case AkR » 1, the cos function oscillates several times 
across the cloud and Ng/Ntot — 1/2 resulting in the loss of contrast. It is worth noting that 
the accumulation of corrections to the wave vectors of the clouds is due to the fact that 
the reflecting pulses do not reverse the clouds' velocities but rather add a constant velocity 
±2^0 to them. This explains the fact that the coherence may be lost due to the presence of 
an external potential even when the nonlinearity is negligible. 

The above-discussed loss of coherence due to incomplete cancelation of the wave vectors 
of the harmonics by the recombination pulse can be also visualized in the following way: the 
wave functions of the V'i harmonics can be represented as'ijj± — ■\/n±{x) exp(i0-|-), where 
is the parabolic phase (the nominal phase ■:L{Mvq/K)x is taken care of by the optical pulses 
and is not included). In the ideal situation, the parabolic phase for each cloud is centered 
at the middle of the cloud. Nonzero values of 5v (or, equivalently, nonzero values of the 
corrections to the wave vectors of the clouds) mean that the phase of each cloud leads or 
lags behind its density envelope. This situation is schematically illustrated in Fig. 1 showing 
the harmonic ■0+ before the recombination with its phase leading the density envelope. We 
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shall show that the optimum recombination corresponds to the situation when the phase 
profiles, not the density envelopes of the clouds are on top of each other immediately before 
the recombination. 

Operation with high values of the contrast can be achieved in several ways. First, the 
relative magnitudes of the nonlinearity and the external potential are adjusted in such 
a way that their effects cancel each other for a given cycle time T (this is not always 
possible). Second, the recombination and/or reflection are conducted with optical pulses 
having different wavelength as compared to the splitting pulse to compensate for the change 
in the wave vectors of the moving clouds. Finally, the recombination is carried out not at 
the nominal recombination time T but at a time such that AkR = 0. The paper is devoted 
to the analysis of the last possibility. 

The rest of the paper is structured as follows: Sec. II provides general formulation of 
the problem. Sec. Ill introduces analj^ical model of the interferometric cycle and Sec. IV 
is devoted to the analysis of the contrast and contains the analytical expressions for the 
optimized recombination time. These expressions are discussed in different limiting cases in 
Sec. V. 

II. FORMULATION OF THE PROBLEM 

The evolution of the condensate in the interferometer in the mean-field limit is described 
by the Gross-Pitaevskii equation 



where '^{r,t) is the wave function of the condensate that is normahzed to one, is the total 
number of atoms, Uq = 47r /i^a^/M characterizes the strength of interatomic interactions, Qg 
is the s-wave scattering length and M is the atomic mass. The potential Vtot = + 
Vopt{t) cos{2kix) is the sum of a confining potential Vsd and an optical potential that is 
created by two counter-propagating laser beams of wavelength X — 27r/ki, which are detuned 
from the atomic resonance to avoid spontaneous emission. The optical potential is used to 
split, rccombine and reverse direction of propagation of the BEG clouds. 
The confining potential is of the form 




(3) 



VMr,t) = V{x,t) + Mu;lrl/2, 



(4) 
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where V{x,t) is slowly-spatially-varying potential due to the environment and is 
the guiding potential providing confinement of the condensate along the two spatial dimen- 
sions r± — {y, z). In the following we shall assume that the condensate is tightly confined 
in the two transverse dimensions and is in the lowest transverse mode of the guide 

Mr±) = ^i^exp(-ri/2ai), (5) 

where a± — {h/Muj±)^/'^ is the transverse oscillator length. Factorizing the wave function of 
the condensate as ^{r,t) — tlj{x,t)t/j±{r±), Eq. (3) can be reduced to the one-dimensional 
equation for the function ip{x,t). Introducing dimensionless coordinate x 2kix and time 
r = t/to, where to — M/^ihkf), this one-dimensional Gross-Pitaevskii equation can be 
written as 

1 



i^ilj{x,T) 



+ v{x, r) -|- ^{t) cos(a;) -|- p\ 



|2 



ip{x,T), (6) 



2dx^ 

where v — {V/h)to, ^ — {Vopt/h)to and p — agN/a^ki. 

The optical potential Q(T)cosa; acts as a diffraction grating for the condensate wave 
function ip. This grating diffracts the condensate into several harmonics separated by mul- 
tiples of the grating wavevector. If the width of the Fourier spectrum of the condensate is 
much smaller than the length of the grating wavevector (one in our dimensional units) , the 
wave function ip{x, t) in Fourier space consists of a series of narrow peaks. It is therefore 
convenient to represent ip{x, t) as 

'^{^: t)^^ i^nix, t) exp(mx) , (7) 

n 

where harmonics' envelopes ipnix) are slowly- varying functions of coordinate as compared 
with the exponentials. The dynamics of these harmonics are governed by the set of coupled 
equations 

n—m+l- (8) 

l,m 

The optical potential ^(t) in Eq. (6) is used to split the initial zero-momentum BEG 
cloud at the beginning of the interferometric cycle into the two harmonics with the mo- 
menta ±1, reverse their direction of propagation in the middle of the cycle and recombine 



them at the end. Dynamics of the BEC due to the optical potential was fully taken into 
account in solving the Gross-Pitaevskii equation numerically. In the analytical model, their 
action was described in terms of simple transformation matrices. For self-consistency of 
the presentation, a brief derivation of parameters of optical pulses used in the numerical 
solution of the Gross-Pitaevskii equation is given in the Appendix. The results presented in 
the Appendix have been previously derived in Refs. [9, 10, 13]. The next section is devoted 
to the development of analytical model describing evolution of the BEC between the optical 
pulses. 



III. PARABOLIC MODEL 



Between the optical pulses the condensate consists of two harmonics with n = ±1 whose 
evolution is described by the set of coupled equations 

^_^d_\^ 1 
dr dx J 2 dx"^ 



(9) 



Introducing the density and phase of each harmonic by the relations 

fp± = v^exp(i0±) 

and using the Thomas-Fermi approximation (neglecting the second derivatives of the density) 
transforms the set of equations (9) to the form 

d ^ d\ - _^ ( ^"f"^ 
dr dx j ^ dx \ ^ dx 



d ^ d 
dr dx 



1 



— V — p{n± + 2n=p). 



2 V dx 

We will describe the external potential v by the first two terms of the Taylor expansion 

1 



(10) 



v{x) = ax -\ — Px 



(11) 



and analyze the set of Eq. (10) in the framework of a parabolic approximation where expres- 
sions for both the density and the pase do not contain terms higher than the second order 
in coordinate: 



3 

8^ 



[X — x±} 



(f)± ^ (p± + K±{X - X±) + -{X - X±)^. 



(12) 



R, x±, K±, (p± and g are functions of time r only. The coefficient 3/8 in tfie expression for 
n± follows from the normalization condition (each harmonics is normalized to 1/2). Note 
that for each cloud its density and phase in Eq. (12) are defined only in the region where 
the density is nonncgativc. Functions x±{t) are positions of the centers of mass of the two 
moving clouds. k± are corrections to their nominal wavevectors (±1) that are due to the 
external potential and the nonlinearity and ip± are the accumulated coordinate-independent 
phases. Finally, R is the half-size of each of the clouds and the parameter g multiplying 
the quadratic part of the phase is analogous to the inverse of the radius of curvature of the 
wavefront of a propagating light beam in optics. 
Using Eq. (12) and the first of Eq. (10), one gets 

i?' = g^R, 

x'^ = ±1 + K±, (13) 

where the prime means differentiation with respect to time. Treatment of the second Eq. (12) 
is slightly complicated by the fact that the regions of existence of n+ and n_ do not coincide. 
Since the functional forms of n and 4> are fixed, the density profile should be projected 
onto n±. To do this, one can choose a set of suitable basis functions defined at the interval 
|C±| < 1) where ^± = (x — x±)/R, that can be used to represent the density n± and the 
phase 0±. The density should then be expressed in terms of the same basis set retaining 
only the functions that describe n± and (j)±. Using Legendre polynomials Pn{C,±) as the basis 
yields 

1 CD 

-^n^{C±)^doTdii±-d2e±, (14) 

where 

do= (^2-^\q\'' + 2\q\'-^\q\^ 0{\q\<2), 
d,^q(^i-3\q\ + ^\qf^ 9{\q\<2), 

d,= (2-^\q\' + 5\qf- l\qA e{\q\ < 2) (15) 



2 8' 

and q = {x+ — X-)/R. The 6'-function in Eq. (15) is equal to one if its argument is a logical 
true and zero if its is a logical false. 

Using Eq. (14) in the second of Eq. (10) yields equations of motion for g, k± and (p±. 
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Combining these with Eq. (13), we get the final set of equations 



R' = gR, 

= ±1 + K,±, 

v'±-Y~ ""^^ ~ ^^"^^ " + ^^^^ 

where prime means differentiation with respect to time. Equations (16) have simple physical 
interpretation. The rates of change of the coordinates of the two clouds are given by 
the relations x'^ = ±1 + k±, i.e., the clouds move with velocities ±1 + k±. The major 
contributions to the velocities ±1 are due to the momenta imparted to the clouds by the 
optical pulses. The corrections k± are due to the external potential (parameters a and (3) 
and the nonlinearity. The cloud "climbing up" the external potential slows down and the one 
moving "downhill" speeds up. If the nonlinearity is repulsive (p > 0), the speeds of the two 
clouds after their separation will be slightly larger than one and if it is attractive, slightly 
less than one. The functions do, di and ^2 given by Eq. (15) describe mutual interaction of 
the two clouds. They depend on the relative displacement of the clouds q — (a;+ — X-)/R 
and are nonzero only when |g| < 2, i.e., when the clouds overlap. The other terms containing 
the nonlinearity parameter p describe self interaction for each of the clouds and are always 
nonzero. 



A. Evolution of k± and x± 

During the interferometric cycle the two BEC clouds il)± may be partially overlapping or 
non-overlapping. In the subsequent analysis, it will be assumed that the size of each cloud 
does not change significantly at the time intervals t oc R that it takes for the clouds to 
pass each other. The conditions of apphcability of this assumption are given by Eq. (29). 
Additionally, it will be assumed that PT"^ <S 1, where T is the duration of the interferometric 
cycle. 
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Time evolution of k± and x± is governed by the set of two coupled equations (cf. (16)) 



3p 



X. 



di, 



.± — ±1 + K,±, 

Solution of Eqs. (17) can be written as 



1 3p 



qo+2T/Ro 



dqdi{q), 



(17) 



(18) 



90 



11 1 

x±(r) = x±,o + (±1 + «±,o)t - -ar^ - -a;±,o/5r^ T g/^r^ 



± 



3p 
"32 



go+2r/ilo 



f'dq'd.iq'), 

J On 



(19) 



where K-i- x±,o and i?o are initial values of x± and R. In deriving Eq. (18) and (19), 
the dynamics of the relative separation between the clouds in evaluating function di was 
approximated by the relation 

?(t) = ?o + ^, (20) 
i.e., the terms with k± were neglected as compared to one in evaluating q{T). 

The interferometric cycle of duration T starts by applying the splitting optical pulses 
to the motionless cloud -00, letting harmonics '0± propagate for the time r/2, reverse their 
directions of propagation by applying the reflection pulses, letting the harmonics %l)± evolve 
for the time T/2 and apply the recombination optical pulses. 

Immediately after the splitting pulses at r = 0, the center of mass of each harmonic 
is x±fl = 0, go = and = 0. The reflection pulse reverses directions of propagation 
of the two harmonics by adding momenta =f2 to the momenta ±1 + k± of ijj±. After the 
reflection pulse the harmonic tp+ becomes ip^ and vice versa. As a result, immediately after 
the reflection pulse, x±fi — x^{T/2) and k,±^o — k^{T/2). 

At the nominal recombination time r — T, the corrections to the velocities k± and the 
center of mass coordinates x± are given by the relations 



K±{T) = -aT±-pT^^p 



1-D,{T/I^) + ^D,{T/Rt) 



(21) 



x^{T)^~aT'±^PT' 



_p 



/ dqD^{q) + —D,{T/Ro) + —D,{T/Rj 

Jt/Rt -"o -n-T 



(22) 
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where 



{2-X+ ^x^) , x<2 



D^{x) = { ^~ " ■ ^ ' (23) 
Di(2) =3/10 , x>2, 

Ro is the size of the harmonics at the separation stage and Rt is the size during recombi- 
nation. 



B. Evolution of g and R 

Evolution of g and R is governed by the set of two coupled equations (see Eq. (16)) 



R' = gR, 

5' = V-/? + ^(l + ^2). (24) 

The explicit expressions for g and R at time intervals r such that R does not change 
significantly, i.e., | Ai?| -C R, are of the form 



g{T) ^g,-pr+-^ dq[l + ^2(9)], (25) 

J go 



R{t) = Ro + Ro 



/3 2 , 3p 
9or--T +^ 



/ dr' / dq[l + d2{q)] 

Jo Jqo 



2 8i?Q JO Jqo 



(26) 



where go — g{0), go — q{0) and Rq — R{0) are initial values of g, q and R. In deriving 
Eqs. (25) and (26), the terms of the order k± in the equations for x± have been neglected 
as compared to one. The dynamics of the relative separation between the clouds in the 
framework of this approximation is given by the expression 

q(r)^qo + ^ (27) 



Equations (25) and (26) are valid provided 



9or, /3t^ ^ « 1 (28) 

In the analysis of Sec.III A and in the rest of the paper it is assumed that the size of each 
cloud does not change significantly during the time t — R that it takes for the clouds to 
pass each other. The conditions of applicability of this approximation are 

goRo, PRl, |- « 1 (29) 
12 



2 

Rn 



Using Eq. (25), we get the following expression for the value of g at the recombination time 
in the limit when R does not change significantly during the interferometric cycle: 



9iT)^go-PT+^^^+^^D,{T/R), 



(30) 



where 



D2{x) = 



tX 



(2 



5™2 _|_ 5™3 
2 ^ \^ 



(31) 



16-^^) ) 2, 

^2(2) = , x>2 

The hmit |Ai?| -C i? can correspond to both r < R when the clouds stay overlapped 
during all the cycle and to t ^ R when they do not overlap most of the cycle. The second 
limit of interest to be considered in this section r ^ R explicitly deals with the situation 
when the clouds do not overlap most of the time. In this limit, the contribution coming from 
the function ^2 in Eq. (24) (interaction between the clouds) can be neglected as compared 
to their self action. The function g in this limit is given by the relation 

1/2 



r dr 



sign go 



r 



9l 



3p 



1 

- - 1 

r 



(3(r' - 1) 



(32) 



where r = R{t)/Rq. Note that Eq. (32) is vahd for any values of R{t). The general 
expression for R{t) can be obtained in terms of elliptic integrals but is too cumbersome to 
be of practical use. In the limit where the relative change in the size of each harmonic is 
small(|r — 1| <^ 1), one gets 

3p 



R{t) ^Ro + Ro 
3p 



gor + 



ARl 



21 



9(^) ^go + 



4RI 



(33) 



These expressions coincide with Eqs. (25) and (26) when t ^ R and describe the situation 
when the clouds do not overlap most of the time but their sizes do not changes significantly 
during all their evolution time. 
In the opposite limit R{t) ^ Rq, 

u 1/2 



g(r) = sign go 



Rq 



9l 



3p 
2^ 



R\t) 
Rl 



(34) 
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C. Evolution of 



In an interferometric experiment, the quantity of interest is not the absolute phase of 
each harmonic </7±, but rather the relative phase Ac/? = — The time evolution of Ac/? 
is governed by the equation 

Ac/?' = - kI) - a{x+ - x_) - \(3{x\ - xl) (35) 

Using results of Section III A and neglecting terms containing products and quadratic or 
higher combinations of a, /3 and p yields 

A(/.(r) = (36) 

IV. THE INTERFERENCE SIGNAL 

The wavefunction of the zero-momentum harmonics after the recombination is given by 
the expression: 

1 1 



tpo{x) = —1= [ip+ix) + ip-{x)] oc —1= \/ n{x — a;+) exp{i9o/2 + iAkx/2) 
v2 v2 L 



+ s/nix + x+) exp(-i6'o/2 - iAkx/2) . (37) 
Here ijj± are the wavefunctions of the ±1 harmonics before the recombination, 

00 = {(p+ - (f-) - {k+x+ - K-X-) + |(a;^ - x'^_) (39) 

and 

Ak^ Ak- gAx, (40) 

where Ak — and Ax — x+ — X-. All quantities in Eq. (39) are evaluated at the 

recombination time. 

If the density envelopes of the ±1 harmonics sufficiently overlap at the recombination 
stage, Eq. (37) can be simplified to 



^o{x) = Vn(x) cos(^o/2 + AA;a;/2) (41) 
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Population of the zero- momentum harmonics A^o is given by the spatial integration of [■^oP 
yielding 

^° ^(1 + ycos^o), (42) 



Ntot 2 

where the contrast of the interference fringes V is given by the expression 

V = r[sin(AA;i?) - AkRcos(AkR)]. (43) 

For AkR <C 1, the population of the zero- momentum state is given by the relation 

No/Ntot^ 008^00/2). 

In this limiting case the population depends on the relative accumulated coordinate- 
independent phase $0 between the two BEC clouds and exhibits interference fringes as 
a function of this phase. 

In the opposite case AkR ^ 1, the cos function in Eq. (41) oscillates several times across 
the cloud and 

No/Ntot = 1/2 

independently of the value of the relative phase shift. 

Equations (22) and (21) show that both the nonlinearity of the condensate p and the 
quadratic contribution to the external potential P can result in nonzero values of Ak given 
by Eq. (40) and thus be responsible for the loss of interferometric contrast as illustrated 
by Fig. 2. This figure shows the contrast V defined by the relation N^/Ntot = (1 + ^)/2, 
where No /Ntot is the relative population of the zero-momentum harmonics at the end of the 
interferometric cycle, as a function of the ratio of the cycle time to the initial size of the 
harmonic T/Rq. The solid line corresponds to the numerical solution of the Gross-Pitaevskii 
equation (6). The dashed line is given by Eq. (43), where AkR is calculated with the help 
of analytical expressions (40), (21), (22) and (30). 

Since the linear slope of the potential is zero {a — 0), 9o — 0. Equation (43) then 
predicts that for AkR — 0, V — 1. As is seen in Fig. 2, the contrast indeed equals one for 
short cycles (small T). Larger values of T correspond to larger interaction times between 
the two clouds and an increase in AkR due to this interaction. As the interaction time 
increases, the contract V given by Eq. (43) goes down from one to small negative value 
resulting in the values of N^/Ntot slightly below 1/2. At times larger than about T/Rq — 1.5 
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the two harmonics completely pass each other and stop overlapping during a part of the 
cycle. The interaction time between the harmonics (the time when they overlap) is now 
smaller than the cycle time and does not depend on it. The contrast and the population of 
the zero-momentum harmonic reach their limiting values. Figure 3 shows the dependence 
of the population of the zeroth-order harmonic after the recombination No/Nfot on the 
relative accumulated phase shift Oq — —aT'^/2. The solid line corresponds to the numerical 
solution of the Gross-Pitaevskii equation (6). The dashed line is Eq. (42) with AkR given 
by Eqs. (21), (22) and (30). As is seen in Fig. 2, the cycle time T/Rq — A corresponds to 
small negative values of the contrast {V ^ —0.2 as given by the Gross-Pitaevskii equation 
and V ~ —0.1 as given by the analytic model). Low values of the contrast result in the 
washout of the interference fringes shown in Fig. 3. It also should be noted that since the 
contrast V is negative, the symmetric recombination with Oq — corresponds not to the 
maximum, but the minimum population A^o of the zero-momentum harmonic. 

Figures 2 and 3 demonstrate that recombination with nonzero value of the linear wavevec- 
tor Ak (see Eq. (40)) washes out the interference fringes. Since Ak is a function of time, 
this effect can be compensated for by conducting the recombination not at the nominal time 
T but at a slightly different time T + AT when AkR = (in general, AT may be both 
positive and negative). Figure 4 shows the contrast V ~ 2No/Ntot — 1, where N^/Ntot is 
relative population of the zero-momentum harmonic, as a function of the time AT. Negative 
(positive) values of AT correspond to the recombination taking place shghtly before (after) 
the nominal recombination time T. The parameters for Fig. 4 are T = 2000, Rq — 500 
and p = 5 with all other parameters being zero. The sohd line is the solution of the Gross- 
Pitaevskii equation and the dashed line is obtained with the help of Eqs. (42), (21), (22) 
and (32). Recombination at the nominal time AT — corresponds to a small value of the 
contrast and a washout of the fringes as is shown in Fig. 3. Figure 4 indicates that if 
the recombination takes place at AT/Rq Ki —0.2, the contrast of the fringes becomes much 
larger. This is confirmed by Fig. 5, which shows No/Nfot versus the relative accumulated 
phase shift = -a[{T + ATf - /2] for A/R = -0.2 and all other parameters the same 
as in Fig. 3. The solid line is the solution of the Gross-Pitaevskii equation and the dashed 
line is the result of the numerical solution of Eqs. (16). 
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The simple estimate using the condition {AkR){T + AT) — yields 



AT 



{AkR)T 
RiAkRy^' 



(44) 



Rt 



The ratio AT/R gives the relative displacement of the two clouds at the recombination time 
T + AT since the clouds pass across each other in time R (each cloud has the size 2R and 
the relative speed is 2). 

The population of the zero-momentum harmonics depends not only on the magnitude of 
AkR, but on the degree of overlap of the two density envelopes at the recombination time 
(see (37)). Estimate (44) takes into account only changes in AkR but not in the overlap 
in evaluating AT. The last can be taken into account in the framework of Eq. (37) at 
the expense of making formulas more cumbersome and turn out to be not very significant. 
As we shall see, Eq. (44) is in a very good qualitative and quantitative agrement with 
the results of numerical solution of the Gross-Pitaevskii equation (6). Finally it should be 
mentioned that nonzero values of AT/R mean incomplete overlap and thus the contrast less 
than one even at the optimized time. The larger is \AT\/ R, the smaller the contrast. The 
estimate (44) implicitly implies that \AT\/R < 1 because correction to the recombination 
time is meaningful only if the clouds overlap at the time T + AT. If, for some set of 
parameters, estimate (44) yields \AT\/R > 1, the coherence can not be recovered for this 
set of parameters. 

Using the explicit expressions for {AkR) and its time derivative obtained with the help 
of Eqs. (40) and (16) results in the relation 



In Equation (45), Ak{T) and Ax{T) are evaluated using Eq. (21) and (22). The function 
g{T) should be evaluated using several different expressions depending on the parameters 
of the problem. For |Ai?| -C -R, g{T) is given by Eq. (30). In this case Rq — Rj' — R. If 
T 3> i?0) 9{T) is given by Eq. (32). Since we are assuming that |Ai?| ^ i? at times it takes 
the clouds to pass through each other, Eqs. (30) and (32) cover all possible situations. If 
\AR\ is not small as compared to i?, the size of the clouds Rt at the end of the cycle should 
be evaluated by numerical integration of Eq. (32). 

Equation (45) is relatively complex because it covers both the case when the size of the 
clouds does not change significantly during the cycle and the opposite limit when the final 



AT Ak - gAx 



(45) 



R 2gR-3pAx/AR^' 
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size is much larger than the initial one. All the relevant physics can be understood by 
discussing the case |Ai?| <S R when Eq. (45) acquires especially simple form 

AT ^ ^ 1 8D,{C)~mVp)e .... 

R 2gR A3C/^ + D2{0+9oRyp-myp)C ^ ^ 

where ( — T/R. 

The contrast at the optimized recombination time can be evaluated by accounting for an 
incomplete overlap of the clouds using Eq. 37 and is given by the approximate expression 



V ^1-- 



3 /atV r, R „ 1 



In — — - + 2 In 2 - - 
AT 2 



(47) 



2 \ R J 

Equations (46) and (47) are the main analytical results of the paper. In Section V, they will 
be analyzed in several illustrative cases. 

V. DISCUSSION 

A. Influence of the nonlinearity p for go = P = 



For l3 = go = 0, Eq. (46) becomes 

AT -2Di{T/R) 



(48) 



R D2{T/R) + 3T/4i? 

Equation (48) shows that the correction to the recombination time depends only on the 
single parameter T/R and does not depend on the nonlinearity of the condensate p (the 
applicability of the parabolic approximation requires pR ^1). This is due to the fact that 
both the A/t and the gR terms are proportional to the nonlinearity parameter p. At small 
values oiT/R, Di{T/R) ^ {3/8){T/Rf, D^iT/R) ^ (3/2)(T/i?) (cf. Eqs. (23) and (31)) 
so 

3 T 

IT = -TOR' 

i.e., AT/R grows linearly with T/R. The correction to the recombination time AT/R 
reaches maximum for T/R ^ 2 when the duration of the cycle is such that the two clouds at 
their maximum separation stop overlapping. At longer cycle times T/R > 2, both Di and 
D2 become constants and AT/ R starts decreasing inversely proportional to T: 

AT 4 /T\"^ , , 
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This behavior has simple physical explanation. The difference between the corrections to 
the propagation velocities of the clouds Ak is due to the nonlinear interaction between the 
clouds and is accumulated only when the clouds overlap (see the definition of Di Eq. (23)). 
For short cycle times T/R < 2, when the clouds overlap during all the cycle, the nonlinear 
efi^ects are accumulated during all times and An oc T^. The parabolic phase described by the 
coefficient g{T) grows linearly with time T, so the correction to the recombination time T 
is a growing function of the cycle time T. For T > R, when the clouds fully separate during 
the cycle, Ak is at its maximum possible value and stops growing further. The quadratic 
phase profile of each cloud, on the other hand, keeps growing as a function of time, i.e., g 
becomes larger, thus resulting in the decrease of AT. 

The dependence of the shift in the recombination time AT/Rq on the cycle time T/Rq 
is shown in Fig. 6 for Rq — 500 and p — 5. The maximum cycle times shown in Fig. 6 
correspond to the maximum separation of the clouds equal to about ten their diameters. The 
dots are the results obtained by direct numerical solution of the Gross-Pitaevskii equation 
(6) and the solid fine is given by Eq. (48). The optimized contrast of the interference fringes 
V at the recombination time T + AT for the parameters of Fig. 7 is shown in Fig. 7. The 
dots correspond to the numerical solution of the Gross-Pitaevskii equation and the solid line 
is given by Eq. (47). For comparison, the dashed line shows the contrast at the nominal 
recombination time T given by Eqs. (43) and (40). The lowest values of the optimized 
contrast V ~ 0.5 correspond to intermediate cycle times T/Rq 2 when the maximum 
separation between the two clouds is equal to their size. Both increasing and decreasing the 
cycle time T improves the contrast. 

Figs. 6 and (7) show that the operation of the atom Michelson interferometer with the 
optimization of the recombination time is possible both in the limit T /2Rq < 1 when the 
clouds overlap during all the cycle and in the opposite limit T/2Rq ^ 1 when the clouds 
are separated most of the time. 

B. Nonzero initial parabolic phase go 

Performing an interferometric cycle with nonzero initial values of the parabolic phase go 
considerably improves the coherence as compared to the case go = provided the sign of go 
is the same as that of the nonlinearity p. The nonzero initial parabolic phase can be acquired 
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by relaxing the confinement frequency lo of the initial trap and letting the condensate evolve 
for some time before the start of the interferometric cycle. Dynamics of the BEC in time- 
dependent parabolic traps in Thomas-Fermi limit has been extensively analyzed (see, e.g. 
[14-16]). In the case of ID expansion corresponding to our situation, the evolution of g{r^ 
and r = R{t)/R{0) is described by the set of equations 

^=^lnr, (51) 

where u!{t) is the trap frequency. The exact value of g depends on the detailed time depen- 
dence of u;{r). Changing u;{t) adiabatically slowly leaves the phase of the condensate flat, 
i.e. g — 0. Since we are interested in the maximum possible value of g, we shall consider 
the limit when the trap frequency is relaxed very fast so that lij{t) = for r > 0. In this 
limit, g is given by the relation (the condensate's initial phase in the trap is zero): 

where Rtr is the initial radius of the condensate in the trap. An extra factor of two in 
Eq. (52) as compared to Eq. (32) is due to the fact that the initial condensate is normalized 
to one whereas the two propagating clouds are normalized to 1/2. 

For the given value of Rtp, g is maximum for r = 3/2. The final size of the condensate 
after the expansion is the initial size Rq of the propagating clouds in the interferometric 
cycle, i.e., Rq/ Rtr — 3/2. The maximum possible value of go is thus given by the relation 

90,max = {^2^ J ' ^^^^ 

In the following we will use the value go — sgo^rnax where the coefficient < s < 1 accounts 
for relaxing the trap with finite speed. 

The correction to the recombination time AT/R given by the equation (46) with P — 
takes the form 

AT ^ -2D,{T/R) 

R {3/A){T/R) + D2{T/R) + s{3R/2py/^ ^ ' 

If the parameter (-R/p) is large, which is typically the case, the corrections to the recombi- 
nation time are small and the contrast is high. This is illustrated by Fig. 8 which shows the 
shift in the recombination time AT/i? as a function of the cycle time T/R using Eq. (54). 

20 



The solid line corresponds to s = when the condensate does not have initial parabolic 
phase. The dashed curve gives AT/R for s — 0.2, when the condensate has been allowed to 
acquire initial parabohc phase. The dots are the results of a numerical solution of the GPE 
with s = 0.2. Figure 8 demonstrates that the shift in the recombination time is considerably 
smaller when the condensate is allowed to expand before the beginning of the cycle. Since 
the two harmonics have larger overlap at the optimal recombination time, the contrast in 
the interference fringes is larger when go ^ 0. 



C. Nonzero parabolic external potential P ^ 

Nonzero values of the parabolic external potential /3 can be due to environment or 
technical imperfections of an experimental apparatus. Equation (46) with (3^0 and go = 
yields 

AT _ 1 8D,{o - myp)e . . 

R 4 3C/4 + D2(C)-(Wp)C' ^ ^ 

where C = T/R. The influence of the parabolic potential on the operation of the atom 

Michelson interferometer is characterized by the parameter b = f3R^/p. Note that since the 

term with f3 in the numerator of Eq. (55) is proportional to the square of the cycle time and 

the denominator grows linearly with time, even small values of b for long enough cycles will 

always result in a complete loss of coherence. 

In the limit of short cycle times C ^ 1; Eq. (54) takes the form 

AT C 3-5 

-R- = -\W^b- ^''^ 

Equation (56) is similar to Eq. (49) but the sign of AT can be both negative and positive 
depending on the value of b. The second difference is in that the coefficient multiplying ( 
may become so large for positive values of 6 5/2, that coherence will be lost even for short 
cycle times. Negative values of P are preferable because they ensure the operation of the 

interferometer at least for short times T/R < 1. If the value of j3 is controlled at the level 
6 <^ 1, the operation of the interferometer is possible for AT/R < 1 and any sign of (3. 
In the limit C > 2, Eq. (55) becomes 

AT_ 12/5 - be , . 

(3 - 46)C ■ ^^^^ 
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If l&l ^ 1, the optimized contrast will be high in the range 




(58) 



If \b\ > 1, the coherence in general will be lost for T/R > 2. 

Additional limitations on the strength of the quadratic potential (3 and the cycle time T 
are due to the fact that the reflection pulses will not operate well if the relative change in 
the velocity of the atomic cloud \Av\/vo exceeds about 0.1. Using the fact that the nominal 
dimensionless velocity Vq ^ 1 and the relation |A^;| = (l/2)|/3|(T/2)2, one gets \(3\T'^ < 0.8. 

D. Recombination at a different wavelength 

The contrast of the interference fringes can be improved by conducting the recombina- 
tion with optical pulses having different wavelength as compared to the splitting pulse to 
compensate for the change in the wave vectors of the moving clouds. The relative change in 
the wavelength of the recombining pulse A A/ A as compared to the separation pulse is given 
by the expression (cf. Eq. (21)) 



As has been discussed in the introduction, the repulsive nonlinearity results in the speeds v of 
the moving harmonics ip^ being smaller than the speed vq imparted by the separation pulse. 
The recombination then should be performed with beams of larger wavelength. Similarly, 
for /3 < (a potential hump) AA > and for /3 > (a potential trough) AA < 0. The 
optimized contrast is determined by the relation 



where Ax is the separation between the centers of the harmonics il)± at the recombination 
time given by Eq. (22). In a typical situation, \Ax\/R -C 1 and the optimized contrast is 
close to one. 
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APPENDIX A: DYNAMICS OF THE BEC DUE TO THE OPTICAL PULSES 



The optical potential is used to split the initial zero-momentum BEC cloud at the begin- 
ning of the interferometric cycle into the two harmonics with the momenta ±1, reverse their 
direction of propagation in the middle of the cycle and recombine them at the end. The 
optical pulses are short and sufficiently intense so that the dynamics of the condensate is 
dominated by the optical potential when the laser beams are on and the diffraction, relative 
displacements of the clouds and the nonlinearity can be neglected. A good quantitative 
description of the BEC dynamics can be obtained keeping only the lowest three harmonics 
with n = 0, ±1 in Eq. (7). The set of Eq. (8) with these approximations reduces to 
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Solution of Eq. (Al) has the form 
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(A3) 
(A4) 
(A5) 
(A6) 



and s — y/1 -\- SVfi. Using Eq. (A2) it is straightforward to show that the momentum reversal 
of the moving BEC clouds %l)±i — > can be achieved with a single pulse of duration 
Tp — 47r and magnitude Q.p — (3/8)^/^. The unitary evolution matrix corresponding to the 
momentum reversal pulse is of the form 



U. 



±1<-*=F1 — 



1 

1 

1 



(A7) 
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Splitting of the zero-momentum cloud ipo into the two harmonics ip±i and the recombination 
(the inverse of the splitting) requires a double pulse sequence. The first pulse with = 
(1/8)^/^ and Tp — 2^/^7r is followed by a period of free evolution when the lasers are turned off 
for a time interval Tev = 2n and then followed by the second optical pulse with = (1/8)^/^ 
and Tp = 2^/^7r. The evolution matrix for the splitting sequence is given by 

-|exp(-i7r/V2) ^ |exp(-i7r/^) 




1 

V2 



(irrelevant common phase has been omitted). 
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V2 



1 
V2 



iexp(— j7r/\/2) 
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FIGURES 
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FIG. 1: The density and phase of the ■0+ harmonic before the recombination. 
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FIG. 2: The contrast V versus the cycle time T/Rq for Rq = 500, p = 5, a = 0, go = and (3 = 0. 
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FIG. 3: Relative population of the zeroth-momentum harmonic Ng/Ntot versus the relative accu- 
mulated phase 9o = —aT^/2 for T = 2000. Other parameters (except nonzero values of a) are as 
in Fig.2. 
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FIG. 5: Relative population of the zero-order harmonic No/Ntot versus the relative accumulated 
phase shift for T/Rq = —0.2. All other parameters are the same as for Fig. 3 
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FIG. 6: The shift in the recombination time AT/Rq as a function of the nominal recombination 
time T/Ro for Rq = 500 and p = 5. The dots are the results of the numerical solution of the GPE 
and the solid line is the analytical model. 



32 




FIG. 7: Optimized contrast of the interference fringes V for the parameters of Fig. 6. The dots are 
the results of the numerical solution of the GPE and the solid line is the analytical model. The 
dashed line is the contrast at the nominal recombination time. 
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FIG. 8: The shift in the recombination time AT/R as a function of the cycle time T/R given by 
Eq. (54) for s = (soHd) and s = 0.2 (dashed curve). The dots are the results of a numerical 
solution of the GPE. For all three curves R/p = 100. 
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